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We investigate slow dynamics of collection of a few antiferromagnetic NiO nanoparticles. The 
zero-field cooled magnetization exhibits size dependent fluctuations. We find memory effects in field 
cooled magnetization, as well as aging effects in thermoremenant magnetization of antiferromagnetic 
nanoparticles. The antiferromagnetic nanoparticles show a stronger memory effect than the corre- 
sponding effect in the ferromagnetic particles, when the distribution of particles include very small 
sizes. The situation reverses for bigger sizes. The relaxation of the magnetization after a sudden 
cooling, heating and removal of fields reiterate the memory effects. We also see a weak signature 
of size-dependent magnetization fluctuations in aging effect of antiferromagnetic nanoparticles. We 
find a two-step relaxation of thermoremenant magnetization in antiferromagnetic case, which differs 
qualitatively from relaxation of ferromagnetic nanoparticles. 
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I. INTRODUCTION 

The interest in the study of magnetism in nanopar- 
ticles has been renewed from the last few decades due 
to their technological^"'^ as well as fundamental research 
aspects4i"— The magnetic properties of nanoparticles are 
dominated by finite-size effects, and the surface anoma- 
lies such as surface anisotropy and roughness.—""— As the 
particle size decreases, the fraction of the spins lying on 
the surface of a nanoparticle increases, thus, making the 
surface play an important role. The reduced coordination 
of the surface spins causes a symmetry lowering locally, 
and leads to a surface anisotropy, that starts dominating 
as the particle size decreases. 

The dynamics of an assembly of nanoparticles at low 
temperatures gained a lot of attention over the last few 
years. In a dilute system of nanoparticles, the inter- 
particle interaction is very small as compared to the 
anisotropy energy of the individual particle. These iso- 
lated particles follow the dynamics in accordance with 
Neel-Brown model^^ and the system is known as super- 
paramagnetic. The giant spin moment of nanoparticles 
thermally fluctuates between their easy directions at high 
temperatures. As the temperature is lowered towards 
a blocking temperature, the relaxation time becomes 
equal to the measuring time and the super spin moments 
freeze along one of their easy directions. As the role of 
interparticle interaction becomes significant, nanoparti- 
cles do not behave like individual particles, rather their 
dynamics is governed by the collective behavior of the 
particles, like in a spin glass i^ri^ This super spin-glass 
phase has been characterized by observations of a crit- 
ical slowing dowuji^""— a divergence in the nonlinear 
susceptibility)^""— and aging and relaxation effects in the 
low-frequency ac susceptibilityjiii The Monte Carlo sim- 
ulations on the system of assembly of nanoparticles show 
aging^ and magnetic relaxation behavior^^ like in a spin 
glass, but simulations of zero field cooling (ZFC) and 
field cooling (FC) susceptibilities show no indication of 
spin-glass ordering)^ The aging and memory effects are 
the two aspects that have been studied extensively in re- 



cent years, however mostly for the case of ferromagnetic 
part icles i^""— i ^ ^ ' ^ ^" "— In this paper our aim is to investi- 
gate these effects for collection of a few antiferromagnetic 
nanoparticles (AFNs) which has received relatively lesser 
attention. 

Neel predicted^i^l that AFNs exhibit weak ferromag- 
netism and superparamagnetism behavior, which is at- 
tributed to a net magnetic moment due to incomplete 
magnetic compensation between the atoms in two sub- 
lattices. These sublattices are identical in every re- 
spect, except that the atomic moments in one sublat- 
tice are antiparallel to that in other. For Nickel Oxide 
(NiO) nanoparticles, the magnetic moment /i predicted 
by Neel's model varies as ~ n3^^j2+, where n is the 
number of spinsi^ Weak ferromagnetism and superpara- 
magnetism were later confirmed by experiment^ on fine 
particles and extremely fine particles of NiO, respectively. 
Mishra and Subrahmanayam'^^ have also concluded a net 
magnetic moment in NiO AFN due to finite size effect in 
Neel-state ordering, which is showing a non-monotonic 
and oscillatory dependence on particle size R. The am- 
plitude of fluctuations were found to be varying linearly 
with i?, consistent with Neel's model. However, the value 
of net magnetic moment due to flnite size effect in Neel- 
state ordering does not quantify the large magnetic mo- 
ment experimentally observed in NiO nanoparticles. This 
discrepancy is due to the ignorance of different order- 
ing of surface spins than the core spins in these studies. 
Monte Carlo simulations on antiferromagnetic particles 
highlighted the dominant role of surface spins in net mag- 
netization of the nanoparticles The surface effects 
have been considered to be the major cause for the large 
magnetic moment in the NiO nanoparticle.— '^i^^'^^''^^ The 
breakdown of the dominant next-nearest neighbor anti- 
ferromagnetic interaction on the surface of the nanoparti- 
cle leads to uncompensated spins. These uncompensated 
spins play a vital role in determining the magnetic be- 
havior of NiO nanoparticles. Thus, an enhancement of 
surface and interface effects make the AFNs an interest- 
ing area of research.— ""—i^""— 

Recently, Mishra and Subrahmanyam^i have shown 
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FIG. 1: A log normal distribution of nanoparticles of sizes 
ranging 1.3ao — 5.3ao- The inset displays the magnetic mo- 
ment vs particle size. Magnetic moment has a non-monotonic 
and oscillatory dependence on R. 



that for NiO nanoparticles the net magnetic moment, 
a combined effect of surface roughness effect and finite- 
size effects in core magnetization, exhibits size depen- 
dent fluctuations in net magnetic moment. These size 
dependent fluctuations in magnetization lead to a dy- 
namics which is qualitatively different from ferromag- 
netic nanoparticles. 

In this paper, we mainly focus on the dynamics of col- 
lection of a few noninteracting NiO nanoparticles using 
master equation approach. We examine the effect of the 
size-dependent magnetization fluctuations on the time- 
dependent properties of sparse assembly of the nanopar- 
ticles. We will compare the ZFC and FC magnetiza- 
tions of antiferromagnetic particles with the ferromag- 
netic case. We will consider the effect of polydispersity 
as well. We will perform a series of heating/cooling pro- 
cesses which were earlier discussed in the case of ferro- 
magnetic nanoparticleSi ^^'^^" — We will also discuss the 
memory effects and aging effects throughly. The organi- 
zation of this paper is as follows. In section [ll] we discuss 
the model. The ZFC and FC magnetizations for various 
distributions are discussed in section lllll and in section lTVl 
we show the memory effects investigations. Aging effects 
has been presented in section [V] Finally we summarize 
in Sec ED 



II. RELAXATION IN SUPERPARAMAGNETS 
AND POLYDISPERSITY 

In an earlier study^^ we have shown that the total 
magnetic moment of NiO antiferromagnetic nanoparticle 
displays size dependent fluctuations as shown in the in- 
set of figure [1] The net magnetic moment shows a trend 
where magnetic moment is very small for smallest size 
particles. Increasing the size, magnetic moment increases 



and reaches a maximum at i? ~ IOoq, and again decreases 
towards the bulk value. On the other hand, the net mag- 
netic moment of ferromagnetic nanoparticles^ shows a 
linear dependence on the size of the particles Hence one 
might expect the role of these size-dependent fluctua- 
tions in magnetization should be manifested in the time 
dependent properties of antiferromagnetic nanoparticles. 

In our simple model, the energy of each particle i is 
contributed by anisotropy energy (either due to the shape 
or the crystalline structure of the particle), and Zeeman 
energy. For the sake of simplicity we assume that the 
direction of field is same as that of anisotropy axes. Thus 



= -Jcv^+n^i^, 



(1) 



where /C is the anisotropy constant and H is the applied 
magnetic field. In the absence of field, the superparam- 
agnetic relaxation time for the thermal activation over 
the energy barrier ICVi is given by t = roexp(/CVi/kBT), 
where tq, the microscopic time is of the order of 10^^ sec 
and fcs is Boltzmann constant. The anisotropy constant 
K. has a typical valued about 4 x 10^^ Jcm^'^. The occu- 
pation probabilities with the magnetic moment parallel 
and antiparallel to the magnetic field direction are de- 
noted by pi (t) and p2 (i) = 1 — pi (t) , respectively. These 
probabilities must satisfy the master equation^ 



A2i(t)(l-PiW), 



(2) 



where the parameters Xi2{t) and A2i(i) are the rate of 
transition of the magnetic moment from the two states 
at time t, given as To"^exp [-/CVi/T(t)] [1 - fiih{t)/T{t)] 
and To-^exp[-/CVi/T(t)][l-FAtih(t)/T(t)] respectively. 



Using the values of parameters Ai2(t) and A2i(<), 
[5] can be simplified as 



m 



Eq. 



(3) 



The magnetic moment of the particle of volume Vi given 
by 

ti{t,V,)^[2pi{t,V,)^l]fx,. (4) 
For mt) = n and T{t) ^ T, we can write 



=M(0,^^z)exp(-t/r) + 



T 



{l-exp(-t/r)}.(5) 



For a constant magnetic field, the above equation gov- 
erns the relaxation of magnetization at each temperature 
step. Thus the total magnetic moment of the system of 
nanoparticles with volume distribution P{Vi) is given by 



fi{t) = / fi{t,V,)P{V^)dV,. 



(6) 



The size-distribution plays a significant role in the over- 
all dynamics of the system of nanoparticles governed by 
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FIG. 2: a)ZFC and FC susceptibilities have been plotted from the solution of Eq. [6] with a heating/cooling rate IO^^tq 
and 10^^ To per temperature step. Ripples can be seen in ZFC susceptibilities. Increasing heating rate lowers the blocking 
temperature, b) ZFC susceptibilities for ferromagnetic and antiferromagnetic cases using same size distribution with heating 
rate lO^^ro has been shown. The susceptibilities are normalized with xfc at T = 0. It can be seen that ripples are absent in 
the ferromagnetic case. 



Eq. [SI As we can see that due to the exponential de- 
pendence of r on the particle size Vi, even a weak poly- 
dispersity may lead to a broad distribution of relaxation 
times, which gives rise to an interesting slow dynamics. 
For a dc measurement, if relaxation time coincides with 
the measurement time scale T„i, we can define^^ a crit- 
ical volume Vb as /CVb = fcB'7Bln(T,n/To), where Tb is 
referred as blocking temperature. The critical volume 
Vb has strong linear dependence on Tb and weakly log- 
arithmic dependence on the observation time scale Tm- 
If the volume of the particle in a polydisperse sys- 
tem is less than Vb, the super spin would have under- 
gone many rotations within the measurement time scale 
with an average magnetic moment zero. These particles 
are termed as superparamagnetic particles. On the other 
hand if Yi > Vb, the super spins can not completely 
rotate within the measurement time window and show 
blocked or frozen behavior. However, the particles having 
volume Vi ~ Vb are in dynamically active regime. The 
systems of magnetic nanoparticles are in general polydis- 
perse. The shape and size of the particles are not well 
known but the particle size distribution is often found 
to be lognormal."*^ We consider the system consisting of 
lognormally distributed, widely dispersed nanoparticles, 
hence non interacting among each other. The volume Vi 
of each particle is obtained from a log normal distribtrtion 



in order to see the effect of size dependent magnetiza- 
tion fluctuation on the relaxation dynamics of assembly 
of nanoparticles. We can solve Eqs. [3l |4] and |6] for any 
heating/cooling process. For example, we can numeri- 
cally solve these equations for a zero field cooled (ZFC) 
protocol. In a genuine ZFC protocol, system is cooled 
from a very high temperature to lowest temperature in 
the absence of magnetic field. Thus system is demagne- 
tized at the lowest temperature. This condition is anal- 
ogous to pi{0) = 1/2 in Eq. ID Now a constant field is 
applied and the system is heated upto high temperature. 
At each temperature change we evolve the system using 
Eq. [5] We can define heating rate in the process as total 
time elapsed at each temperature change. Thus heat- 
ing/cooling rate IO^^tq per temperatue unit corresponds 
to heating/cooling process in which system is relaxed for 
t — IO^^tq at each temperature step. 

In this paper, volume Vi is measured in the units of 
average volume V. Also the average anisotropic energy 
ICV is taken as the unit of energy, which is equal to 480K 
for V = 193ao'^. Thus by setting fcs = 1, we use ICV as 
a unit of temperature T and and field H. Hereinafter we 
use a dimensionless quantity h = /isH/A^flQ as a unit of 
field, e.g., h = 0.01 is equivalent to a magnetic field 300 
Gauss. 



P(V,;a;u) = 



1 



ctVV^tt 



exp 



-(In(Vi) 



2(72 



(7) 



III. ZFC AND FC MAGNETIZATIONS 



where v — ln{V), V is the mean size and a the width 
of the distribution. The distribution consists of Iff* par- 
ticles of sizes between R = I.Sgq and R = 5.3ao, where 
ao(= 4.I7A) is the lattice parameter of NiOJ^ The to- 
tal number of particles are purposely chosen to be small 



The relaxation phenomenon of nanoparticles is often 
investigated using two different protocols, ZFC magneti- 
zation and FC magnetization measurements. In a ZFC 
magnetization measurement, the system is first demagne- 
tized at a very high temperature and then cooled down to 
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FIG. 3: ZFC susceptibility for various size distributions is 
shown. The role of size dependent magnetization fluctuations 
can be explicitly seen as a ripple in the curve. For bigger sizes 
distribution R = 9.3ao — IS.Sao, the curve is smoother. 



little role of magnetization fluctuation as comapared to 
smaller sizes, which indeed is the case for the size range 
(9ao — 13ao) as shown in Fig. [31 We find that for the size 
range (9ao — 13ao), ZFC curve is almost smooth, while 
for lower sizes, it exhibits ripples. 

During an FC measurement, the system is cooled in 
the presence of a probing field from higher temperatures 
to a low temperature. We find that the FC susceptibil- 
ity Xfc coincides with xzfc at higher temperature but 
departs from ZFC curve at lower temperatures, however 
well above the blocking temperature, and tends to a con- 
stant value with further lowering the temperature. The 
blocking temperature shows a substantial dependence on 
the heating rate. For infinitely slow heating rate, Tb 
approaches to zero and the ZFC curve shows similar be- 
havior as FC curve. We also find that FC magnetization 
never decreases as the temperature is lowered which is a 
characteristic feature of superparamagnetsi^ 



IV. MEMORY EFFECT 



a low temperature in a zero magnetic field. A small mag- 
netic field is then applied and the magnetization is calcu- 
lated as a function of increasing temperature. We have 
shown a plot of FC-ZFC magnetization with temperature 
for polydisperse NiO nanoparticles in Fig. [UJa) using 
Eqs. [31 m and m with heating/cooling rates 2.4 x IO^^tq 
and lO^^To per temperature step. For a ZFC process we 
find that increasing the temperature susceptibility xzfc 
increases, attains a maximum value at a blocking tem- 
perature (Tb), and then starts decreasing. We see that 
the average blocking temperature depends on the heat- 
ing rate. Increasing the heating rate lowers the blocking 
temperature. For example, for 2.4 x IO^^tq per temper- 
ature unit heating rate, blocking temperature is 0.053, 
whereas for heating rate IO^^tq per temperature unit, 
blocking temperature is 0.041. 

In Fig. [2{b), we have shown the comparison of ZFC 
magnetization for antiferromagnetic particles with the 
ferromagnetic particles of same sizes, where susceptibility 
has been normalized by FC value at T = 0. The over- 
all ZFC magnetization behavior is same for both cases 
except the presence of ripples in ZFC magnetization for 
antiferromagnetic particles. These ripples are attributed 
due to the size dependent fluctuations in magnetization. 
As we increase the heating rate these ripples become 
more pronounced. This can be seen in Fig. [5^, where 
increasing the heating rate enhances the ripple in magne- 
tization. We have also investigated the effect of polidis- 
persity by incorporating various size distributions. We 
consider three distributions dist.A, dist.B and dist.C in 
which sizes are in the range of 1.3ao — 5.3ao, 5.3ao — 9.3ao 
and 9.3ao — 13.3ao respectively. The mean sizes in these 
distributions are 3.5ao, 7.5ao and ll.Sao respectively. We 
can argue from Fig. [31 that the smaller size distribu- 
tions exhibit more effects of fluctuations in magnetiza- 
tions than the bigger sizes. For bigger sizes, we expect a 



Recently Sun et al}^ have reported a striking result 
showing memory effects in the dc magnetization by a se- 
ries of measurements on a permalloy NigiFeig nanopar- 
ticle sample. These measurements include FC and ZFC 
relaxations under the influence of temperature and mag- 
netic field. They cooled the sample in 50 Oe field at 
a constant cooling rate of 2 K per minute from 200 K 
to Tbaso — 10 K. After reaching at Tbase, the sample 
was heated continuously at the same rate upto 200 K. 
The M{T) curves thus obtained are the normal FC curve 
which is referred by them as reference curve. The sam- 
ple is cooled again at the same rate, but with temporary 
stops at T = 70 K, 50 K and 30 K below blocking tem- 
perature Tb with a wait time tw = 4 hours at each stop. 
During each stop, the field was also turned off to let the 
magnetization relax. After each pause the magnetic field 
was reapplied and cooling was resumed. The cooling pro- 
cedure produces a steplike M (T) curve. After reaching 
Tbasc, the sample is warmed continuously at the same rate 
to Th in the presence of the 50 Oe field. The M(T) curve 
thus obtained also shows the steplike behavior around 
each stops. Sun et al suggested that this 'memory ef- 
fect' indicates the possibility of hierarchical organization 
of metastable states resulting from interparticle interac- 
tions. Since hierarchical organization requires a large 
number of degree of freedom to be coupled, the mem- 
ory effect may not arise due to the thermal relaxation of 
independent particle. However more recently, Sasaki et 
al?- and Tsoi et al^ have reported similar results for the 
noninteracting or weakly interacting superparamagnetic 
system of ferritin nanoparticles and Fe203 nanoparticles 
respectively. In these studies, the dynamics of the system 
of nanoparticles are assumed to be governed by a broad 
distribution of particle relaxation times arising from the 
distribution of particle sizes and sample inhomogeneities. 
Experiments on NiO nanoparticles by Bisht and Rajee\^ 
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FIG. 4: Memory effect in the dc magnetization with a stop at T = 0.039. Memory effects shown in (a), (b) and (c) correspond 
to antiferromagnetic nanoparticles for dist.A (1.3ao — 5.3ao), dist.B (5.3ao — 9.3ao) and dist.C (9.3ao — 13.3ao) with a stop at 
T = 0.039, while (d), (e) and (f) correspond to ferromagnetic nanoparticles of respective size distributions. 



also confirms a weak memory effect in these particles. 

Chakraverty et a^- have investigated the effect of poly- 
dispersity and interactions among the particles in an as- 
sembly of nickel ferrite nanoparticles embedded in a host 
non magnetic Si02 matrix. They found that either tun- 
ing the interparticle interaction or tailoring the particle 
size distribution in nanosized magnetic system leads to 
important application in memory devices. 

We perform a similar study in thermoremenant- 
magnetization (TRM) protocol as that by Sun et ali^, as 
shown in Fig. |31Ja)-(h). Figures Slja), (b) and (c) corre- 
spond to system of antiferromagnetic nanoparticles with 
distributions dist.A, dist.B and dist.C and FigsSJ^d), (e) 
and (f) correspond to ferromagnetic nanoparticles of size 
distributions dist.A, dist.B and dist.C. We first cool the 
system from a very high temperature to Tbasc = 0.0018 
with a field h ~ 0.01 and then again heat to get the 
Ref. curves in FigsU] (a)-(h) with a cooling/heating rate 
2.4 X lO^^To per temperature step. We again cool the sys- 
tem from a high temperature to Tbase but with a stop of 
lO^ro at T = 0.039. The field is cut during the stop. Af- 
ter the pause the field is again applied and the system is 
again cooled up to base temperature Tbase- The process 
is shown as T J, curves in Figs|Ua)-(h). Finally we heat 



the system at the same rate as that of cooling without 
any stop, shown as T f curves in Figs |31[a)-(h). We find 
that the magnetization shows an upturn exactly around 
T — 0.039. Further heating recovers the Ref. curve. 

We can also shed some light on the role of poly- 
dispersity on memory effect. We compare the relative 
strength of memory effect in these distributions by in- 
troducing a parameter, memory fraction i.e., the ratio 
of AM/A/^'=f at the stop during memory measurement, 
where AM — M — AI^"^. In all the cases, we fix the wait- 
ing time at stop to be IO^'^tq. The calculated values of 
the memory fraction for antiferromagnetic case is 13.3%, 
36.1% and 27.39% respectively. For ferromagnetic case 
the memory fraction takes values 11.4%, 18% and 25% 
respectively. An interesting observation in the case of an- 
tiferromagnetic nanoparticle is the higher value of mem- 
ory fraction for dist.B than that for dist.A and dist.C. 
As we know that the finite-size and surface roughness 
effects are responsible for the enhancement of net mag- 
netization for the intermediate size ranges dist.B, the net 
memory dip is more in this case. However, for ferromag- 
netic nanoparticles, where the net magnetic moment en- 
hances linearly with the size, we find an increasing trend 
of memory fraction. One more interesting aspect is that 
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FIG. 5: (a) ZFC relaxation and (e) FC relaxation for a system of antiferromagnetic nanoparticles. Relaxation curves during ti 
and on lograthimic time scale for negative-temperature, field-change and positive-temperature cycles for ZFC case are shown 
in (b), (c) and (d) respectively and for FC case are plotted in (f), (g) and (h) respectively. 



the memory fraction of antiferromagnetic nanoparticles 
are more than that of ferromagnetic nanoparticles in all 
the cases. Thus using the scale of memory fraction, we 
conclude that for smaller size distributions, memory ef- 
fect is stronger in antiferromagnetic nanoparticles than 
ferromagnetic particles. However, it is obvious that for 
larger sizes distributions, memory effect in ferromagnetic 
particles exceeds to their antiferromagnetic counterpart. 

We have examined the magnetization relaxation in 
NiO nanoparticles by a series of heating and cooling pro- 
cesses to understand the memory effects as discussed by 
sun et alr^ These relaxation studies are performed under 
both the ZFC and the TRM protocol. 

In a ZFC protocol for the distribution dist.A, as shown 
in Fig. [5ja), the sample is cooled down to Tq — 0.036 
in a zero field and a magnetic field h = 0.01 is ap- 
plied after zero wait time. The magnetization relax- 
ation is calculated as a function of time for time inter- 
val ti = lO^^'ro . After ti interval, we employ following 
three distinct routes, either changing temperature or ap- 



plying magnetic field to study the relaxation for time 
t2 = 1 X 1014 - 4 x 10"to. 

A - Negative-temperature cycle: the sample is 
quenched in constant field to a lower temperature, 
Tiow = 0.024, and the magnetization is recorded for 
t2 time interval. 

B - Field-change cycle: the sample is quenched to 
a lower temperature, Tiow — 0.024 with cutting 
off the magnetic field, and the magnetization is 
recorded for t2 time period. 

C - Positive-temperature cycle: the sample is 
quenched in constant field to a higher temperature, 
Tiiigh — 0.048, and the magnetization is recorded 
again for t2 time period. 

At last we bring the system back to initial tempera- 
ture To = 0.036 and a constant field h — 0.01 and find 
the relaxation for interval — 2 x 10^'' — 3 x IO^^'tq. 



7 



We find that during negative-temperature cycle for in- 
terval t2, the particles which were dynamically active at 
To becomes frozen at lower temperature Tiow = 0.036 
and the smaller particles which should be dynamically 
active at Tiow have already been polarized during inter- 
val ti , hence the relaxation becomes very weak during t2 
and a flat relaxation can be be observed in the flgure. 
When the temperature returns to Tg, the magnetization 
also comes back to the level it reached before the tem- 
porary cooling. Moreover, by plotting the data points 
during and ti in Fig. [SJb), we find a continuity be- 
tween two regions. During the field-change cycle, we find 
that the relaxation during t2 is fast with opposite sign 
but as we return back to Tq with a field h = 0.01 again 
applied, the magnetization again comes back to the level 
as before temporary cooling and field change and the 
relaxation curves during ti and are continuous (Fig. 
[SJc)). These results validate the memory effects and also 
show that the relaxation at lower temperatures has no 
influence on the states at higher temperatures. Finally 
during positive-temperature cycle, raising the tempera- 
ture enhances the number of dynamically active particles 
which are taking part in relaxation process. The parti- 
cles which were frozen at Tq become dynamically active 
at increased temperature and polarize themselves during 
t2 interval. As we return back to initial temperature To, 
these particles remain frozen in polarized state. Hence 
the relaxation curves during ti and ^3 show an abrupt 
jump in magnetization as shown in Fig. [SJd). 

We have also performed the relaxation studies in TRM 
protocol. Here, as shown in Fig. EJe), the sample is 
cooled down to Tq = 0.036 in a magnetic fleld h = 0.01. 
The fleld is cut off for a zero wait time and the mag- 
netization relaxation is calculated as a function of time. 
After a time interval ti = lO^^ro, again we employ fol- 
lowing three different processes to study the relaxation 
for t2 = lO^Vo. 

A - Negative-temperature cycle: the sample is 
quenched in zero field to a lower temperature, 
Tiow — 0.024, and the magnetization is considered 
for interval t2 = ti = IO^^tq. 

B - Field-change cycle: the sample is quenched to a 
lower temperature, Tiow = 0.024 with the magnetic 
field h = 0.01 applied, and the magnetization is 
considered again for a time ^2 = ^1 = IO^^tq. 

C - Positive-temperature cycle: the sample is 
quenched in zero field to a higher temperature. 
Thigh = 0.048, and the magnetization is considered 
for a time ^2=^1 = IO^^'tq. 

At last we quench all the above three processes sepa- 
rately to To = 0.036 with zero field to observe the relax- 
ation for interval 2 x lO^'' — 3 x IO^^tq. We have also 
shown the magnetization during ti and for the three 
different processes occurred during ^2 in Fig. [Iff), E^g) 
and \5lh). We find that the logarithmic relaxations in 



Fig. [5jf) and Fig. ^g) are continuous which again in- 
dicates the memory effect. On the other hand during 
the positive-temperature cycle, as shown in Fig. EJh), 
the magnetization is not continuous. This is because the 
bigger size particles, which were frozen at To, become 
dynamically active at higher temperature Thigh. These 
particles depolarize themselves during t2 interval and re- 
main frozen in depolarized state, when temperature is 
brought back to lower temperature To. 

The above relaxation studies in ZFC and TRM pro- 
tocols reveal that the negative-temperature cycle and 
field-change cycle show memory effect, while positive- 
temperature cycle does not imprint memory during the 
cycle. Thus the memory effect is only due to the fast re- 
laxation of smaller particles which respond to temporary 
cooling and field change. 

V. AGING EFFECT 

Aging effect is a well studied phenomenon in spin glass 
system.™!^^ Recently it has gained a lot of attention in 
the system of nanoparticles, where slow dynamics study 
becomes important to characterize both the superspin 
glass behavior and the superparamagnetism .^'^^i^°'^^i^^ 
Most of the experiments are performed by measurement 
of time-dependent ZFC and TRM magnetizations. In a 
TRM protocol the system is cooled in a field to a base 
temperature Tbasc below blocking temperature Tb. After 
a waiting time ^wi the magnetic field is switched off and 
one observes the relaxation in magnetization. It is found 
that the time dependence of the magnetization depends 
on the waiting time t^. These studies also show that even 
a noninteracitng system of ferromagnetic nanoparticles 
can exhibit the aging effect though weak. 

We have studied the aging effect using a polydisperse 
system of NiO nanoparticles in TRM protocol. Our in- 
vestigation is carried out by cooling the system in the 
presence of magnetic field h = 0.01 upto base tempera- 
ture Tbase = 0.024 and cutting the field off after a wait 
time <w to let the system relax. For all the three size 
distributions, system of nanoparticles, either ferromag- 
netic or antiferromagnetic, show wait time dependence. 
In Figsinja), (b) and (c), we plot the magnetization of an- 
tiferromagnetic nanoparticles with logarithmic time scale 
for different wait times IO^^'to, IO^^to and IO^^tq for dis- 
tributions dist.A, dist.B and dist.C respectively. We 
also plot relaxation curve for the case of ferromagnetic 
nanoparticles using same size distributions dist.A, dist.B 
and dist.C in FigstB^d), (e) and(f). We see that relax- 
ation of antiferromagnetic nanoparticles is qualitatively 
different than that of ferromagnetic particles. The time 
dependence of thermoremenant magnetization in ferro- 
magnetic nanoparticles can be described by a stretched 
exponential function 

/i(t)=A*(0)exp(-(t/T„)"), (8) 

where r™ is response time which has a correspon- 
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FIG. 6: Aging effect for waiting times 10^ to, 10 to and 10^ to. (a), (b) and (c) correspond to antiferromagnetic nanoparticles 
for dtst.A (1.3ao — 5.3ao), dist.B (5.3ao — 9. Sao) and dist.C (9. Sao — 13. Sao), (d), (e) and (f) correspond to ferromagnetic case 
for respective size distributions. The solid curves in ferromagnetic cases are fits to Eq. |8]and those in antiferromagnetic cases 
are fits to Eq. [5] 



dence with the peak position of relaxation rate S{t) — 
9Af/91ogjg(t) versus logjg(t) curve. 

Best fit of Eq. [8] to aging data of distribution dist.A 
shows Tm — 1.078 X IO^^tq and n = 0.0566 for wait- 
ing time lO^-^To; r„ = 1.105 x IO^tq and n = 0.0569 
for waiting time IO^^tq; and Tm — 1.2053 x IO^^'tq and 
n = 0.0582 for waiting time IO^^tq. Using same equa- 
tion to fit the aging data of distribution dist.B reveal 
Tm = 1.001 X IO^^tq and n = 0.115 for waiting time 
lO^ro; Tm = 1.025 x 10"ro and n = 0.116 for wait- 
ing time lO^^To; and Tm = 1.12 x 10"ro and n = 0.119 
for waiting time IO^^tq. For distribution dist.C con- 
sisting of bigger sizes, the fitting parameters are t„i — 
1.0 X 10"to and n = 0.192 for waiting time IO^tq; 
Tm — 1.028 X IO^'^tq and n = 0.194 for waiting time 
lO^Vo; and r„ = 1.128 x 10"ro and n = 0.199 for wait- 
ing time IO^^tq. 

For all the distributions discussed above, we find a 
small increase in parameter n with waiting time t^. Thus 
parameter n can be useful to quantitatively describe ag- 
ing effect. 



However, Eq. [5] does not satisfy the relaxation of an- 
tiferromagnetic nanoparticles. As we can see multiple 
peaks in relaxation curves shown in Figs. [Ha), (b) and 
(c), one step exponential decay can not be sufficient to 
describe the behavior in this case. We define a two step 
stretched exponential decay for magnetization of the form 

^,it) = Ai(O) [exp(-(t/r^J"0 +exp(-(t/r„,J"^)] . (9) 

The best fit of Eq. [5] to magnetization vs logj^Q (t) data 
for various size distributions are listed below. For dist.A, 
we find the parameters Tmi = 1.03 x IO^'^tq, ni = 0.0124, 
Tm2 = 7.282 X lO^Vo and n2 = 0.392, for waiting 
time 10"ro; r^^ = 1.03 x lO^Vo, ni = 0.0157, r^^ = 
7.40x lO^^To and n2 = 0.385, for waiting time IO^^tq; and 
Tm, = 1.05 X lO^To, ni = 0.0158, r^^ = 6.52 x IO^Vq 
and rt2 — 0.373, for waiting time IO^^tq. A slight increase 
in ni and a decrease in 712 with wait time can be seen for 
this size distribution. 

For dist.B, fitting parameters are Tm, = 3.48 x IO^^tq, 
m = 1.603, Tm^ = 1.12 X lO^Vo and n2 = 0.183, for 
waiting time IO^^tq; Tm, = 3.46 x IO^'^tq, ni = 1.604, 



9 




\og^M\) log,i^Ao) logio(^/'^o) 



FIG. 7: Relaxation rates for waiting times 10^*ro, lO^^ro and W^^'tq. (a), (b) and (c) correspond to antiferromagnetic 
nanoparticles for dist.A (1.3ao — 5.3ao), dist.B (5.3ao — 9. Sao) and dist.C (9.3ao — 13.3ao). (d), (e) and (f) correspond to 
ferromagnetic case for respective size distributions. A peak can be seen in the ferromagnetic cases, which depends on wait 
time tw and corresponds to Tm in Eq. [8] Muhiple peaks can be seen in antiferromagnetic cases and the relaxation of TRM 
magnetization is governed by a two step stretched exponential function given by Eq. 



Tm^ = 1-0 X IO^^tq and n2 — 0.181, for waiting time 
lO^Vo; and r^^ = 3.53 x IQ^^tq, ni = 1.749, = 
1.0 X IQi^To and na = 0.193, for waiting time IO^^tq. We 
see that for dist.B, waiting parameters rti and n2 show 
increasing trend with waiting time. 

However for bigger size distribution dist.C, fitting pa- 
rameters are = 1-01 x IO^'^tq, ni — 0.34, t„i^ — 
1.25 X lO^^To and n2 — 0.619, for waiting time IO^^tq; 
Tm^ = 1.02 X 10"ro, ni = 0.35, r^^ = 1.22 x IO^^to 
and 71-2 = 0.613, for waiting time IO^^tq; and 7:,„^ — 
1.3 X lO^To, ni = 0.329, r™, = 1.19 x lO^^ro and 
71-2 = 0.606, for waiting time IO^^tq. We find that pa- 
rameters Til and 712 are decreasing with increasing wait 
time. 

This may easily be understood from a weak depen- 
dence of blocking volume on logarithmic observation 
time. For the case of ferromagnetic nanoparticle, de- 
magnetization of particles with logarithmic time is linear 
but the same is not true for antiferromagnetic case. We 
also find that magnetization decays faster for bigger size 
distribution dist.C than smaller sizes in both case ferro- 



magnetic as well as antiferromagnetic nanoparticles. 



VI. CONCLUSIONS 

We have studied the effect of size-dependent magneti- 
zation fluctuations on the dynamics of the polydisperse 
system of AFNs by solving two state model analytically. 
A collection of a few antiferromagnetic nanoparticles has 
been has been Numerical calculation of ZFC magnetiza- 
tion shows ripples in the curve which is absent in ferro- 
magnetic particles of same size. These ripples are sig- 
nature of size-dependent fluctuations in magnetization 
and they become more pronounced as heating rate is in- 
creased. The distribution of sizes also play an important 
role in the time dependent properties of the polydisperse 
system of nanoparticles. Ripples in ZFC magnetization 
curve are more highlighted for smaller size distribution 
and disappear for larger sizes. A broad distribution of 
particle relaxation times arising from the polydispersity 
is found to be responsible for the dynamics of the system 
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of nanoparticles viz superparamagnetic. The memory ef- 
fect and a weak aging effect has also been observed in a 
noninteracting polydisperse assembly of nanoparticles for 
various size distributions. For very small nanoparticles, 
Memory effect is more in antiferromagnetic case than fer- 
romagnetic case. The situation reverses for bigger size 
nanoparticles. We have also discussed various relaxation 
measurements with sudden cooling, heating and removal 
of fields to validate the memory effects. We have found 
that the striking memory effects in system of AFNs are 
indeed originated from polydispersity. A fitting to aging 
data in antiferromagnetic nanoparticles shows a two step 
stretched exponential decay, as contrast to ferromagnetic 
case, where magnetization show a stretched exponential 
decay. This can also be confirmed by observed multiple 
peaks in the relaxation rate versus log2o(t) curve. In the 
case of ferromagnetic nanoparticles, aging parameter n 
shows an increasing trend with waiting time. For anti- 
ferromagnetic nanoparticles, depending upon the distri- 



bution of sizes, aging parameters ni and n2 can increase 
or decrease with waiting time. 
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